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Summary 

The content of this Quarterly Report can be summarized as follows: 

The analysis of the rotational dynamics of the satellite has been focused 
on the rotational amplitude increase of the satellite, with respect to the 
tether, during retrieval. The dependence of the rotational amplitude upon the 
tether tension variation to the power 1/4 has been thoroughly investigated. The 
damping of rotational oscillations achievable by reel control has also been 
quantified while an alternative solution that makes use of a lever arm attached 
with a universal joint to the satellite has been proposed. Comparison simula- 
tions between the SAO and the Martin Marietta (MMA) computer code of retrieval 
maneuvers have also been carried out. The agreement between the two, completely 
independent, codes has been extremely close, demonstrating the reliability of 
the models. 

The slack tether dynamics during reel jams has been analytically investi- 
gated in order to identify the limits of applicability of the SLACK3 computer 
code to this particular case. Test runs with SLACK3 have also been carried out. 


^RECEDING PAGE NOT 


1 -^ 



Page 4 


Eigure 

Figure 


E igure 
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Figure 


Captions 


1. Wobble of the subsatellite with no initial angular velocity 
and an angular displacement halfway between the x and y 
principal axes. Part (a) Motion of the tip of the body z- 
axis as seen from the wire for the first 25 seconds. Part 

(b) Motion of the tip of the body z-axis for the first 90 
seconds. Part (c) The rotation angle of the body x-axis 
relative to the inertial y-z plane. Part (d) The component 
of angular velocity along the body z-axis. Part (e) The 
kinetic energy vs. time. 

2. The relationship between the inertial axes and the body axes. 

3. A 600 sec test run for retrieval of the subsatellite with the 
same initial conditions for rotation as Figure 1. Part (a) 
Tether length vs. time. Part (b) Tension vs. time. Part 

(c) In-plane angle vs. time. Part (d) Out-of-plane displace- 
ment of the tether attachment point vs. time. Part (e) 

Angular velocity component along the body x-axis. Part (f) 

Angular velocity component along the body y-axis. Part (g) 

Angular velocity component along the body z-axis. Part (h) 

The motion of the body z-axis as seen from the wire. Part 

(i) The angle 6 between the wire and the body z-axis. Part 

(j) The rotation angle of the subsatellite about the direction 
of the wire. Part (k) The rotational kinetic energy vs. time. 

4. Retrieval simulation for 10,000 seconds of orbital time with a 
wobble of the subsatellite. Part (a) Tether length vs. time. 
Part (b) Tension vs. time. Part (c) In-plane libration angle. 
Part (d) Angu-lar velocity component along the body z-axis. 

Part (e) Rotation angle of the subsatellite about the body 
z-axis. Part (f) The angle between the wire and the radius 
vector to the attachment point. Part (g) The rotational 
kinetic energy vs. time. 

5. Simulation of z-axis spin control using program ROTAT with 
stiffness and damping reduced by a factor of ten in the 
control algorithm. Part (a) Component of the angular 
velocity along the body z-axis. Part (b) Rotation angle 
of the subsatellite about the z-axis. Part (c) The work 
done by the z-axis attitude control thruster. Part (d) 

The integrated torque of the z-axis attitude control 
thruster . 
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F igure 


F igure 
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Figure 


F igure 
Figure 


Figure 


Simulation of a single axis rotation of the subsatellite 
using program ROTAT with an exponential tension decrease. 
Part (a) The tension vs. time. Part (b) The angle 8 
between the wire and the vector to the attachment point vs. 
time. Part (c) The rotational kinetic energy vs. time. 

Part (d) The subsatellite rotation angle in inertial 
coordinates vs. time. 

7. Simulation of a circular wobble using program ROTAT with 
an exponential tension decrease. Part (a) Tension vs. 
time. Part (b) The angle 8 between the wire and the 
radius vector to the attachment point vs. time. Part 
(c) The rotational kinetic energy vs. time. Part (d) 

The motion of the tip of the body z-axis as seen from 
the direction of the wire. 

8. Simulation of a single axis rotation with program ROTAT 
with a sinusoidal tension variation. Part (a) The tension 
vs. time. Part (b) The angle 8 between the wire and the 
radius vector to the attachment point vs. time. Part (c) 
The rotational kinetic energy vs. time. Part (d) The 
subsatellite rotation angle in inertial coordinates vs. 
time. Part (e) The angular velocity vs. time. 

9. Simulation of a single axis rotation using program ROTAT 
with a sinusoidal variation of the tether angle. Part (a) 
Tether angle vs. time. Part (b) Rotation angle of the 
subsatellite in inertial space. Part (c) Rotation angle 
with respect to the tether. Part (d) Component of the 
angular velocity along the body x-axis. Part (e) The 
rotational kinetic energy vs. time. 

10. Rotation angle vs. time for the first 565 seconds of 
retrieval starting at 20 km. 

11. Kinetic energy vs. time for rotation of the subsatellite 
with wire damping included in the model . 

12. Lever arm for damping rotational oscillations of the 
subsatellite. Part (a) Lever arm remains aligned with 
the wire. Part (b) Lever arm does not follow the 
direction of the wire closely. 

13. Retrieval from 20 km with a single axis rotation of the 
subsatellite and two terms in the tension control law. 

Part (a) Tether length vs. time. Part (b) In-plane 
tether angle vs. time. Part (c) Tension vs. time. 

Part (d) Rotation angle of the subsatellite with respect 
to the tether vs . time . 
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Figure 14. Retrieval from 20 km without rotational dynamics and three 
terms in the tension control law. Part (a) Tether length 
vs. time. Part (b) In-plane angle of the tether vs. 
time. Part (c) Tension vs; time. 

Figure 15. Retrieval from 20 km with a single axis rotation of the 
subsatellite and three terms in the tension control 
law. Part (a) Tether length vs. time. Part (b) 

In-plane angle vs. time. Part (c) Tension vs. time. 

Part (d) Rotation angle with respect to the tether vs. 
time for 10,000 seconds. Part (e) Rotation angle vs. time 
for the first 500 seconds. Part (f) Rotation angle vs. time 
for the last 500 seconds. Part (g) The rotational kinetic 
energy vs. time. 

Figure 16. Snapshot of the strain and velocity profiles following a 
reel jam. 

Figure 17. Final tether velocity V f as a function of reel velocity 
at the moment of j am , V . 

Figure 18. SLACK3 simulation of reel jam. Retrieval angle 45° , 

1 km deployed tether at jam. 
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1.0 INTRODUCTION 

This is the seventh Quarterly Report submitted by SAO under contract NAS8- 
36160, "The Investigation of Tethered Satellite System Dynamics," Dr. Enrico C. 
Lorenzini , PI and covers the period from 15 February 1986 through 14 May 1986. 

• 

2.0 TECHNICAL ACTIVITY DURING REPORTING PERIOD AND PROGRAM STATUS 

2.1 Subsatellite Rotational Dynamics 

2.1.1 Initial Conditions Eor Retrieval With A Wobble Of The Subsatellite - 

Section 2.2 of Quarterly Report #5 shows simulations of the buildup of 
angular oscillations during retrieval for a one dimensional rotation with the 
rotational angular velocity parallel to the orbital angular velocity. The simu- 
lations showed that the rotational energy at the end of retrieval was only about 
1/6 of the initial energy. The angular amplitude increased by more than a 
factor of 2. Another retrieval simulation has been done with unequal moments of 
inertia and a two-dimensional wobble. The wobble alternates between circular 
and linear motion because of the asymmetry. No exact equations have been de- 
rived to calculate initial conditions for a circular wobble with unequal moments 
of inertia. However, the fact that the motion goes through a linear phase 
provides an alternate method of setting up initial conditions. If the subsatel- 
lite is given an initial rotation about an axis that is half way between x and y 
axes, this generates an oscillation with both x and y components. Since the 
frequencies of the components are different the oscillation alternates between 
linear and circular phases. 
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A short test run has been done to demonstrate that a circular wobble can be 
created from an initial displacement having both x and y components. The x and 
y axes are taken to be the axes of the principal moments of inertia. A rotation 
about only x or y should not produce a circular wobble. The initial conditions 
are expressed as Euler angles <p , 8, and ip which are rotations about the body x, 
z, and x axes respectively. The tether is connected to the tip of the body z- 
axis . The initial conditions <p = 0 , 8 = 20°, and ip = 45° should produce the 
desired result (the value of <p is immaterial since the 8 rotation is about the 
body x-axis) . A run has been done for 25 seconds with output every .3 seconds. 

Figure la shows the motion of the tip of the body z-axis as seen from the 
wire. The coordinate axes on the graphs are the inertial y and z axes with the 
wire along the -x axis. The relationship of the axes is shown in Figure 2. The 
motion starts at the left side of the Figure la. The motion has become nearly 
circular by the end of the run. The run has been extended to 90 seconds and 
Figure lb shows the wobble pattern over that time. The plot shows an extremely 
regular pattern covering a square area. Figure lc shows the rotation angle of 
the body x-axis relative to the inertial y-z plane. The motion is small because 
of the small amplitude of the wobble. Figure Id shows the component of the 
angular velocity along the body z-axis. Figure le shows the kinetic energy vs. 
time. The initial value is zero since the initial conditions are for a dis- 
placement with no angular velocity. The motion is initially linear with a large 
exchange between kinetic and potential energy. As the motion circularizes at 
around 30 seconds the variations decrease. The small variation in rotation 
angle seen in Figure lc is an indication of the low coupling between the wobble 
about the x and y axes and the spin about the z-axis. The body could have a 
spin about the z-axis in addition to the slight rotation back and forth produced 


by the wobble. 
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The initial conditions from the test run with program ROTAT have been put 
into DUMBEL for a short retrieval test run. The initial conditions for the 
attitude dynamics are 9 = 2°, <f> = 0 , ip = 45° and 9 = <j> = ip ~ 0 . The initial 
conditions for the tether have no initial in-plane or out-of-plane libration and 
no retrieval velocity. A run has been done for 60 0 seconds with output every 
second. Eigure 3 gives the results. Parts a, b, c, and d give the tether 
length (cm) , tension (dynes) , in-plane tether angle (deg) , and out-of-plane 
tether displacement (cm) . The large in-plane tether displacement masks the 
variations due to the wobble. The effects of the wobble show up clearly in the 
out-of-plane tether displacement. Parts e, f, and g give the components of the 
angular velocity along the body x,y, and z axes. The period of the oscillations 
about the x axis is shorter because of the smaller moment of inertia. The 
period depends on the tension which is decreasing during the runs. Parts h, i, 
and j give the motion of the z-axis as seen from the wire, the angle between the 
z-axis and the wire (9) , and the rotation angle about the wire direction. The 
z-axis plot is somewhat noisy because the wobble period is on the order of 10 
seconds and the output interval is 1 second. The plot of 9 shows the periods of 
linear and circular oscillation. The average value of 9 is increasing during 
the run. The back and forth motion of the spin angle has a small amplitude but 
a systematic effect is developing. Eigure 3k shows the kinetic energy. The 
dominant feature is the alternation between linear and circular oscillation. 
There is a slight decrease in the average kinetic energy during the run. No 
thrusters are used during the run. The initial kinetic energy of rotation is 
589.9 ergs. This is due solely to the orbital angular velocity. The body has 
no initial angular velocity with respect to the tether coordinate system. Sub- 
tracting the work done by the wire from the kinetic energy gives a value which 
is constant to within a range of about 589.6 to 589.9 ergs. The changes in 
kinetic energy are solely due to the work done by the torque of the wire. 
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Figure 3a+ 


Figure 3bi 
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Figure 3ci Figure 3di 
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2.1.2 Subsatellite Orientation With Respect To The Tether Coordinate System - 

There are four coordinate systems that are useful in describing the dynam- 
ics of the tethered subsatellite. The equations of motion are integrated in an 
inertial coordinate system since they have their simplest form in this system. 
The position and velocity of the Shuttle can be used to define an orbital coor- 
dinate system which is useful in describing the libration angles of the tether. 
The principle axes of the subsatellite form a coordinate system. The relation- 
ship of the body axes to the inertial axes are specified by 9 direction cosines 
which are integrated vs. time to determine the rotational dynamics. The orien- 
tation of the subsatellite is influenced primarily by the torque exerted by the 
tether. It is useful therefore to describe the orientation of the subsatellite 
with respect to the tether coordinate system. 

The tether coordinate system can be defined in terms of the velocity v of 
the Shuttle and the force F exerted by the tether on the subsatellite. The z- 
axis of the tether system is taken to be in the direction of the wire force F. 
That is 


z = F/| F | (1) 


The x-axis is taken to be in the direction of the vector X given by 

X = z x v (2) 


The unit vector for the x-axis is 


x = X/| X 


( 3 ) 



The unit vector in the direction of the y-axis is 


y = z x x 
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(4) 


A vector P in inertial coordinates can be transformed to the tether coordinate 
system using the equation 


P" = B P 


(5) 


where P" is the vector in the tether system. In the body axis system the vector 
is P ' given by 


P' = A P 


Premultiplying equation (5) by B T gives 


B T P" = P 


Substituting equation (7) into equation (6) gives 


( 6 ) 


(7) 


P* = AB T P" 


( 8 ) 


Defining 


C = A B T 


( 9 ) 


we have 
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P ' = C P" (10) 

The matrix C defines the orientation of the body with respect to the tether 
coordinate system. The Euler angles <p , 6 and ip can be computed from the ele- 

ments of the matrix C. The matrix B T is 


u 

Yi 

Zi^ 


Y2 

Zj 

l X 3 

y 3 

Z 3 J 


( 11 ) 


where Xi , y A , and Zi are the components of the unit vectors x, y, and z respec- 
tively. The validity of equation (11) is obvious from equation (7) . For exam- 
ple if 


P" 


(l\ 

0 

w 


Then 



Similar relations hold for y and z. 


( 12 ) 


(13) 


A subroutine called ROTANG has been written using the methods developed in 
equations (1-4) (9) , and (11) . This subroutine has as input the state vectors 
of the Shuttle and subsatellite and the matrix A of direction cosines. The 
output is the matrix C and the Euler angles <p , 9, and ip. This subroutine has 
been added to program RSTAVEC which processes the file of state vectors created 
by the numerical integration programs such as DUMBEL. Figures 3h , 3i , and 3j 
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were produced using the output of subroutine ROTANG. 

The format of the output for program DUMBEL has been changed to eliminate 
various quantities which were not generally useful and add some new output. All 
rotational dynamics quantities given in the inertial coordinate system have been 
eliminated. At each output point the Euler angles of the subsatellite are 
printed with respect to both the orbital coordinate system and the tether coor- 
dinate system (using the new subroutine ANGROT) . File FOROll has been changed 
to contain the Euler angle 6 in the tether system rather than in the orbital 
system. As a test, the case of Figure 3 has been rerun for 250 seconds with the 
new version of DUMBEL. The results agree with the previous runs of DUMBEL and 
postprocessing by RSTAVEC. 


2.1.3 Retrieval With A Wobble Of The Subsatellite - 

The run of Figure 3 has been done for 10,000 seconds of orbital time with 
output at 1.0 second intervals. The Shuttle of mass 100 metric tons is at 296 
km with a 500 kg subsatellite deployed upward on a 20 km, 2.54 mm diameter 
tether. The initial Euler angles of the subsatellite in the tether coordinate 
system are 8 = 2°, <p = 0 , and ip = 45°. The satellite radius is 80 cm and the 
moments of inertia are 80, 96, and 99 kg-m J for the x, y, and z axes respec- 
tively. The tether is initially aligned with the local vertical. The tension 
is regulated by the retrieval tension control algorithm according to equation 
(2.1.16) of Quarterly Report #3. 

Figure 4 shows the results of the retrieval simulation. Parts (a), (b) , 
and (c) give the tether length, tension, and in-plane libration angle. The 
control law damps out the in-plane librations during the retrieval . Parts (d) 
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and (e) show the z component of angular velocity and the spin angle about the z- 
axis . The wobble introduces a small spin but the spin rate is very small and 
the satellite does not do a full revolution even without, any spin control. 
Figure 4f shows the angle between the wire and the radius vector to the attach- 
ment point. The amplitude of the wobble goes from 2° in the beginning to about 
4.8° at 10,000 seconds. The motion alternates between linear and circular os- 
cillation because of beats between the x and y components of the wobble. Figure 
4g shows the kinetic energy as a function of time. The maximum kinetic energy 
goes from about 2.7 x 10 s ergs at the beginning to about 4.8 x 10 4 ergs at the 
end. It is difficult to determine exact values from the plot because the 1.0 
second output interval does not give very good resolution. 

The period of the wobble lengthens as the tension decreases. The frequency 
w of the wobble is, for small amplitudes, 

w = i/ rF/I (14) 

The initial value of the tension is 5.35 x 10 6 dynes. With I x = 80 kg-m 2 and I y 
= 96 kg-m 2 , and r = 80 cm, equation (14) gives a period of 8.59 seconds for the 
x-component. The average observed periods are about 8.56 and 9.40 seconds re- 
spectively. The motion is complex with unequal moments of inertia so that the 
periods of the motions are not exactly constant. The average period of the 
wobble computed with I = 88 kg-m 2 is about 9 seconds. Table 1 lists the ten- 
sion, and average computed wobble period at 1000 second intervals during the 


run . 
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Fi8ure 4at Figure 4b* 
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Time 

Tension 

Wobble 

(sec) 

(dynes x 10 6 ) 

Period (sec 

0 

5.35 

9.0 

1000 

3.80 

10.7 

2000 

2.35 

13.6 

3000 

1.48 

17.1 

4000 

1.15 

19.5 

5000 

.79 

23.5 

6000 

.53 

28.6 

7000 

. 38 

33.7 

8000 

. 27 

40 .3 

9000 

.18 

48.6 

10 ,000 

. 13 

57.8 


The numerical integration of the retrieval run took about one hour of com- 
puter time. Processing of the output files took anywhere from 5 minutes for 
shorter files to 20 minutes for the state vector file processed by RSTAVEC. 


Since the spin of the subsatellite induced by the wobble is very small it 
should be possible to maintain the orientation of the subsatellite using the z- 
axis attitude thruster. A simulation has been done using program ROTAT with 
both a restoring term and a damping term in the thruster control algorithm. The 
restoring stiffness is 1/10 that of the wire and the damping coefficient is 1/10 
of the value for critical damping of a wobble under the restoring torque of the 
wire. The initial conditions are 6 — 2°, <p - .6025040423 rad/sec and ip = 
-.6021370130 rad/sec. Figure 5 shows the results of a 25 second run with output 
every .3 seconds. Parts (a) and (b) are the z-component of the angular velocity 
and the spin angle about the z-axis. The maximum angular excursion is about .16 
degrees. Parts (c) and (d) show the work done by the thruster, and the inte- 
grated torque of the thruster. The work done is positive because the thruster 
is being used in a spring mode, but the amount of fuel used is small because the 
wobble has a low coupling with spin about the z-axis for a 2° wobble amplitude. 

















INTEGRATED TORQUE (DYNE-CM-SEC) Z-AXIS THRUSTER WORK (ERGS) 
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The decrease in rotational energy during retrieval shown in Figure 4g is 
due solely to the work done by the wire. Subtracting the work done by the wire 
from the rotational kinetic energy gives a constant value. The only damping in 
the model is the damping term in the tension control law used for retrieval. 
The damping term Fp in the algorithm is 

E d = 4mn (£-£ c ) (15) 

where m is the subsatellite mass (500 kg) , fl is the orbital angular velocity 
(.001158 rad/sec) , t is the reeling velocity and i c is the "commanded" reeling 
velocity. For the given values of m and Q, the retrieval damping coefficient b r 
is 


b r s= 4 m fl = 2316 dynes/ (cm/sec) (16) 

If the retrieval is following the commanded profile, F D is zero. A wobble of 
the subsatellite with a 2° amplitude produces a maximum wire stretching d£ of 

d£ = r (1-cos 2°) = .049 cm (17) 

This occurs in a quarter of a cycle, or about 2.5 seconds for a 10 second 

period. The average velocity is .02 cm/sec. The damping force is therefore 
about 45 dynes and the power about .9 ergs/second. The rotational kinetic en- 
ergy decreased by about 2.2 x 10 5 ergs in 10,000 seconds, or an average of 22 

ergs/second. From this crude analysis it appears that the damping in the con- 

trol law is more than an order of magnitude too small to account for the loss of 
rotational energy observed during the retrieval run. 
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The energy integrals show clearly that the loss of rotational energy during 
retrieval is due solely to the work done by the wire. In general, the wire can 
do either positive or negative work on the subsatellite depending on how the 
tension varies with time. Any asymmetry of the tension during one wobble cycle 
can alter the rotational energy. The principle term in the retrieval algorithm 
is proportional to the length of the wire and should not produce the type of 
asymmetry that would alter the energy. 

The next section describes a mechanism whereby the rotational energy can 
change as a result of slowly varying changes in the tension. 


2.1.4 Angular Oscillation Amplitude As A Function Of Tether Tension - 

Suppose a spherically symmetrical subsatellite is executing a perfect wob- 
ble such that the attachment point of the wire moves in a circle perpendicular 
to the mean direction of the wire. If the tension slowly decreases, the radius 
of rotation of the attachment point will increase. In the process, the body 
does work against the restoring torque of the wire. Let us assume that the 
tension changes are slow enough that the attachment point continues to rotate in 
a circle whose radius varies as the tension changes . Under these assumptions it 
is possible to derive an analytic expression for the angle of the wobble as a 
function of the tether tension using the principle that the decrease in rota- 
tional kinetic energy is equal to the work done against the torque of the wire. 

The final equation obtained from the derivation is 


/v 


f-l/4 


(18) 
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where 9 is the angular amplitude of circular wobble and F is the wire tension. 
The complete details of the derivation are contained in Section 1.0 of Technical. 
Note TP86-002, January 1986. This technical note is included as Appendix B of 
Quarterly Report #6, March 1986 for this contract. 

Equation (18) has been used to analyze the case of Figure 4. Table 2 shows 

as a function of time the wire tension F, the wobble amplitude 9, and the 

predicted wobble angle 9 a (F 0 /F) 1/4 . The last column is the percentage error in 

the observed wobble amplitude. The percentage error is increasing with time, 

/ 

but in general the data agrees quite well with equation (18) . The discrepancies 
observed could be due to the effect of the damping in the retrieval control law. 
As a check on this, the rotational kinetic energy E of the last entry has been 
computed using the formula 


E = r F 9 3 / 2 (19) 

where r is the radius of the subsatellite. Using this equation, the actual 
energy of the last entry is 38,100 ergs and the predicted value is 41,300 ergs. 
The discrepancy is about 32 0 0 ergs which is .3 ergs/sec. This is well within 
the estimate of .9 ergs/sec for the damping in the retrieval law. The damping 


should be greatest at the beginning when the frequency is highest. 
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Table 2 


t 

(sec) 

F 

(dynes x 10 6 ) 

e 

(deg) 

0 O (F O A) 1/4 

%Error 

0 

5.356 

2 . 000 

2 . 000 

. 0 

101 

4.817 

2.044 

2 . 054 

.5 

208 

4.434 

2.096 

2.097 

. 0 

485 

4.000 

2 . 145 

2.151 

.3 

1008 

3.795 

2.172 

2.180 

.4 

2045 

2 . 288 

2.463 

2.474 

.4 

3074 

1 . 445 

2.760 

2.775 

.5 

4086 

1.119 

2.931 

2.958 

.9 

5025 

.780 

3.194 

3.238 

1.4 

6020 

.526 

3.520 

3.573 

1.5 

7047 

. 376 

3.808 

3.885 

2 . 0 

8054 

. 261 

4.153 

4.255 

2.5 

9033 

.182 

4.520 

4.658 

3.1 

9912 

.134 

4.831 

5.029 

4.1 


Equation (18) was derived under very specialized conditions. In the re- 
trieval run of Figure 4, there are various factors that differ from the condi- 
tions under which the equation was derived. The damping in the control law is 
probably the principal factor causing the observed discrepancies. The retrieval 
included both linear and circular wobble of the subsatellite as a result of the 
unequal moments of inertia. The relative success of the formula suggests that 
it is valid for a linear as well as a circular wobble of the subsatellite. 


As noted earlier the long retrieval run used to create the data of Figure 4 
required about a hour or computer time. In addition, some of the postprocessing 
runs took as much as 20 minutes. Almost all of the information needed for 
plotting is available during the integration. In order to save time, it has 
been decided to add extra output files to the integration program which contain 
various data ready for plotting. A file FOR012 has been added containing the 
time, Euler angle 6 in the tether coordinate system, the rotational kinetic 
energy E k , and the "total energy" obtained by subtracting the work done by the 
wire from the kinetic energy. Another file FOR013 has been created containing 
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the time, and the components of the body x and z axes in the tether coordinate 
system. This file is used to plot the motion of the axes, and can be used to 
compute the spin angle of the subsatellite. A new program called PLOTMO has 
been created for generating printer page plots of any of the quantities on the 
new output files. 

SI 

Since the retrieval run involved various effects that differ from the con- 
ditions under which equation (18) was derived it was decided to run some tests 
under more controlled conditions using program ROTAT. In the first of these 
tests the tension E instead of being computed from the stretch of the wire is 
varied according to the equation 


E = E 0 e~ ct (20) 

where F 0 is the initial tension and c is the time constant. A test run has been 

done with the parameters 9 = 2° , <j> = ip = Q , and c = .01. The integration was 

done for 200 seconds. Figure 6 shows the results of the run. Part (a) is the 
tension vs. time, part (b) is the Euler angle 6, part (c) is the rotational 

kinetic energy and part (d) is the spin angle of the subsatellite. The Euler 

angle in part (b) is the absolute value of the spin angle in part (d) for this 
case. The tension in part (a) shows more noise than expected. The problem 
turned out to be an error in the call to the subroutine that computes the values 

of various parameters at each output time. The subroutine was given the inte- 
gration time instead of the output time. The tension in this case is an ex- 
plicit function of the time rather than being a function of the state vector 

(which was correctly interpolated and passed to the subroutine) . 

Table 3 shows a comparison of the observed oscillation amplitudes with 


equation (18) . 
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Table 3 


t 

(sec) 

F 

(dynes x 10 7 ) 

9 

(deg) 

0 O (F O /F)V4 

/(Error 

0.000000 

. 3990530 

2.000000 

2 .000000 


15.506974 

. 3417312 

2.079027 

2.079057 

. 0014 

32.317990 

. 2888520 

2.168278 

2.168297 

. 0009 

50 .673044 

. 2404143 

2.270112 

2.270113 

. 0001 

70.884693 

.1964185 

2 . 387769 

2.387768 

- . 0001 

93.371668 

.1568638 

2.525848 

2.525847 

- .0001 

118.711907 

.1217507 

2 . 691027 

2.691038 

. 0004 

147.736041 

.0910796 

2.893544 

2.893559 

. 0005 

181.701447 

.0648502 

3.149947 

3.149995 

. 0015 


The time in the first column is the time of a maximum in the plot of the angle 9 
or the spin angle. These times are obtained by quadratic interpolation of the 
output points at .3 second intervals. The wire tension E in the second column 
is calculated from the time t using equation (20) . The values of 9 in the third 
column are computed by quadratic interpolation from the output at .3 second 
intervals. The fourth column is the angle 8 computed from equation (18). The 
last column is the percentage error in the computed value of 8 . The agreement 
of the observed values of 9 with equation (18) is virtually exact in this simu- 
lation. The observed residuals may be interpolation error. The close agreement 
indicates that equation (18) can be applied to the case of a linear oscillation. 

A second test has been run applying equation (20) to a circular wobble. 
The initial conditions are 9 = 2°, ^ = .5680462585788 rad/sec and = 
-.5677002201461 rad/sec and the moments of inertia are 99 kg-m J for all axes in 
order to give a perfect circular wobble initially. The run was done with and 
without the error in the time passed to the output subroutine. The output value 
of F was the only quantity affected by the error. Figure 7 shows the results of 
the run. Part (a) is the tension vs. time. The plot is free of the slight 
noise observed in Figure 6a. Part (b) is the angle 9 vs. time, part (c) is the 
rotational kinetic energy, and part (d) is the motion of the body z'-axis as 
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seen from the wire. 


The plot of the angle 8 shows an oscillation superimposed on the secular 
change. In order to determine the phases of the oscillation, program PLOTMO has 
been modified to compute and plot the derivative as well as the value of each 
quantity being displayed. This facility has been used to locate the maxima and 
minima of the slope of the 9 plot. If we assume the fluctuation in 9 is sinus- 
oidal, a maximum or minimum of 8 should correspond to a mean value of 9 . The 
point 1/4 cycle after a maximum of 6 should be a maximum value of 9 and a point 
1/4 cycle after a minimum of 9 should be a minimum of 9 . The data of Figure 7 
covers a 200 second interval with output every .3 seconds. A set of points near 
the end at various phases of the 9 oscillation has been compared to equation 
(18) . The results are shown in Table 4 


Table 4 


t 

Phase 

F 

9 

M F„/F)V« 

%Error 

(sec) 

of 9 

(dynes x 10 6 ) 

(deg) 

(deg) 


0 . 0 


3.9905303 

2.00000 



167.4 

AV 

.7482068 

3.03923 

3.03936 

.004 

170.7 

MAX 

.7239189 

3.07142 

3.06454 

- .224 

174.0 

AV 

.7004195 

3 . 09000 

3.08993 

- .002 

177 . 3 

MIN 

.6776828 

3.10881 

3.11552 

.216 

180.6 

AV 

.6556843 

3.14119 

3.14133 

.004 

183.9 

MAX 

.6343998 

3.17444 

3.16736 

- .223 

187.5 

AV 

.6119676 

3.19660 

3.19599 

- .019 


The results show that the average value of 9 agrees well with equation (18) . 
The peak values of 8 are higher than the calculated value and the minimum values 
are lower. Table 5 shows a comparison with equation (18) of various average 
values of 9 selected from the extrema of 9 . 
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Table 5 


t 

F 

6 

9 0 (F 0 /F) V4 

/(Error 

(sec) 

(dynes x 10 6 ) 

(deg) 



0 . 0 

3.9905303 

2 .00000 



20.4 

3.2541273 

2.10480 

2.10465 

- .007 

33.0 

2.8688869 

2.17195 

2.17200 

.002 

46.5 

2.5065921 

2.24653 

2.24655 

.001 

64.8 

2.0874101 

2 . 35164 

2.35172 

. 003 

80.7 

1.7805532 

2.44709 

2.44708 

. 000 

102.6 

1.4303571 

2.58505 

2.58480 

- . 010 

143 . 4 

.9511572 

2.86241 

2.86236 

- . 002 

174.0 

.7004195 

3.09000 

3.08993 

- . 002 

: agreement 

is quite good over 

' the whole 

range. The 

errors observed 

■ to the lack of interpolation 

to obtain 

the values at 

the extrema of 


times can be in error by up to .05 seconds. 


In the next test of equation (18) , the tension F is varied according to the 
equation 


F = F 0 [1 — .9 sin (.Olmt)] 


( 21 ) 


For convenience in plotting, two extra quantities have been added to the output 
file FOR012. They are the wire tension and wire length making a total of 6 
quantities in all. The initial conditions for the run are 0 = 2°, <f> = ip = 0, so 
that the oscillation is linear with no wobble. Figure 8 shows the results of 
the 200 second simulation. Part (a) is the tension vs. time, part (b) is the 
Euler angle 9, part (c) is the rotational kinetic energy, part (d) is the rota- 
tion angle, and part (e) is the angular velocity. The angle 9 in part (b) is 
the absolute value of the rotation angle in part (d) . The angular velocity 
about the body y and z axes is zero since the rotation is one dimensional . 
Table 6 shows a comparison of the peak values of 9 with equation (18) . The 
first column is the time in seconds, the second is the tension calculated from 
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equation (21) using the time in the first column, the third column is the peak 
value of 0, the fourth column is the value of B calculated from equation (18), 
and the last column is the percentage error in the calculated value. 

Table 6 


t 

F 


8 6 

0 (F 0 /F)V4 

%Error 


(sec) 

(dynes x 10 6 ) 


(deg) 

(deg) 



0 .00000 

3.99053 


2.00000 




28.47810 

1.18918 


2.70097 

2.70692 

.22 


57.00310 

.48562 


3.41381 

3.38620 

- .81 


70 . 01141 

1.08572 


2.78049 

2.76922 

- .41 


99.19274 

3.89946 


2 . 02477 

2 . 01158 

- .65 


123.46641 

6.40483 


1.78884 

1.77689 

- .67 


148.27146 

7.57671 


1.71533 

1.70380 

- .67 


177.14195 

6.35358 


1.79243 

1.78046 

- .67 


196.14038 

4.42494 


1.96189 

1.94899 

- .66 


The tension in 

Figure 8a has 

its 

lowest value at 

50 seconds . 

The period of 

the 

oscillations is 

longest at this 

time . . Equation 

(18) was derived under the 

as- 

sumption that the tension varies 

slowly over one 

cycle. This 

condition is prob- 

ably not well 

satisfied in 

the 

interval around 

50 seconds. 

In Table 6 

the 


greatest percentage errors occur in the range around 50 seconds. From 100 sec- 
onds on, the percentage error remains nearly constant. The period of the oscil- 
lations is shortest in the range 100 to 200 seconds because the tension is 
largest. The amplitude of the oscillation varies according to equation (18), 
becoming larger as the tension decreases, and becoming smaller as the tension 
increases . 

A factor that can affect the rotational energy of the subsatellite is 
changes in angle of the tether. This occurs during retrieval as a result of 
Coriolis forces on the subsatellite. A test run has been done using program 
ROTAT in order to see the effect on the oscillation amplitude and rotational 
kinetic energy of the subsatellite as a result of changes in the angle of the 
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tether. In this run the angle a of the tether has been varied according to the 
equation 


a = .2 sin (.01 n t) radians 


( 22 ) 


The amplitude of .2 radians is equal to 11.46 degrees which is on the order of 
the retrieval angle. The run has been done for 200 seconds with output every .3 
seconds. The initial conditions are 8 = 2°, <f> = \p — 0. Figure 9 gives the 
results of the run. Part (a) is the tether angle vs. time, part (b) is the spin 
angle of the subsatellite in inertial space, part (c) is the spin angle with 
respect to the tether, part (d) is the angular velocity along the body x'-axis, 
and part (e) is the rotational kinetic energy. The energy in part (e) is pro- 
portional to the square of the angular velocity in part (d) . The angle of the 
subsatellite follows the angle of the tether. The amplitude of the oscillation 
with respect to the tether in part (c) is nearly constant. Table 7 gives some 
peak values of the rotational angle 9 at around 0, 50, 100, 150, and 200 sec- 
onds . 


Table 7 


t 

(sec) 

.5407 
6.0765 
44.8041 
50 . 3374 
55.8702 
94.5978 
100.1341 
105.6633 
144.3947 
149.9275 
155.4597 
188.6554 
194.1911 


8 

(deg) 

2 . 0981 
-2.1053 
2 . 0639 
-2.1338 
2 . 0641 
-2.1046 
2.0988 
-2 . 0924 
2.1333 
-2 . 0635 
2.1333 
2.1109 
-2.0923 
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Figure 9at 


Figure 9b 4- 
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The average value of the peak amplitudes is 2.0995 deg ± .0251 deg. The change 
in angle of the tether clearly has an effect on the amplitude of the rotation. 
The rate of change of the tether angle is much greater in this test case than 
the rates during retrieval . 

Equation (18) has been used to analyze the data from two 10,000 second 
retrieval runs described in Quarterly Report #5, December 1985. These runs are 
described beginning at the bottom of page 15 and the results are plotted in 
Figures 7 and 8. It was noted at that time that there was a significant energy 
loss which was assumed due to the damping in the control law. Figure 7 of the 
referenced report shows results of a retrieval run starting with equilibrium 
initial conditions for retrieval at a tether length of 2 km. The oscillation of 
the subsatellite is one-dimensional with no coning motion. The initial angle of 
rotation of the subsatellite is 8 deg and the initial angle of the tether is 

12.1846 deg so that the angle of the subsatellite with respect to the tether is 

4.1846 deg. Since the subsatellite oscillates about the angle of the tether it 
should oscillate between about 8 deg and 16.3692 deg as measured in the orbiting 
reference frame. Table 8 of the present report shows a comparison of the oscil- 
lation angle with equation (18) . The equation agrees very well with the actual 
oscillation angles. The values of the oscillation amplitude are obtained by 
quadratic interpolation to find the peaks and the times given in the first 
column are the times of the peaks. The output values of the tension F are 
printed at 3 second intervals and the values listed are a rough interpolation. 
The lack of accurate interpolation probably contributes to the errors listed. 
There does not appear to be any systematic loss of amplitude such as observed in 
Table 2. The rotation frequency is much lower for the case of Table 8 since the 
run was started at 2 km. This reduces the effect of the damping term in the 


control algorithm. 
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Table 8 


t 

F 

e 

MFo/F) 1 /* 

/(Error 

(sec) 

(dynes x 10 s ) 

(deg) 

(deg) 


42.2 

4.042 

4.2229 



993.0 

2.853 

4.6047 

4.6071 

. 05 

2019.4 

1.985 

5.0499 

5.0446 

- .10 

2987.3 

1.401 

5.5020 

5.0357 

.03 

3987.5 

.9804 

6.0153 

6.0173 

. 03 

5000.2 

.6850 

6.5786 

6.5816 

.05 

5989.9 

.4820 

7.1786 

7.1861 

.10 

7044.0 

.3325 

7.8779 

7.8851 

. 09 

7986.7 

.2391 

8.5528 

8.5627 

.11 

8911.4 

.1737 

9.2686 

9.2748 

. 07 

9774.9 

.1293 

9.9799 

9.9851 

.05 


Figure 8 of Quarterly Report #5 gives the results of a retrieval run 
started at 20 km. Part (d) of the Figure shows the amplitude of the rotation of 
the subsatellite as measured in the tether coordinate system. This amplitude 
has been compared to the amplitude predicted by equation (18). as a function of 
the tension. The comparison showed errors of 5.8, 14.1, 13.3 and 13.7 percent 
at 1005, 2502, 4995, and 9941 seconds, respectively. The damping in the control 
law is not large enough to account for the results. The discrepancies appear to 
be due primarily to the bug described on page 34 of the same report. Due to an 
error in variable names, the torque applied to the subsatellite contained only 
the l and not the ( l~l c ) term in the control law. Since the retrieval was not 
started in equilibrium, the (i-l c ) term is significant. In the 2 km retrieval 
the term is negligible since the run uses equilibrium initial conditions 
for retrieval. Since the erroneous torque is proportional to £, the amplitudes 
of Figure 8d have been compared to equation (18) with F replaced by t. This 
gives percentage errors of 1.2, 2.5, 3.7, and 4.1 at 1005, 2502, 4995, and 9941 
seconds respectively, which are more reasonable. 
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A short rerun has been done for the case of Figure 8 in Quarterly Report #5 
using the current version of the computer program. Table 9 shows a comparison 
of .the oscillation amplitude with equation (18) . The results agree very well 
with the formula. Eor the original run the percentage errors were 2.0, 4.0, 
5.2, 5.9 and 6.0 at 98.3, 203.5, 296.5, 403.8 and 498.9 seconds respectively 
because of the mistake in the torque calculation. Figure 10 shows a plot of the 
rotation angle vs. time during the 565 second rerun. The angle is increasing 
because of the rapid decrease in tension as the subsatellite accelerates. 


Table 9 


t 

F 

6 

*o(F 0 /F)V4 

%Error 

(sec) 

(dynes x 10 6 ) 

(deg) 

(deg) 


16.36 

5.78 

2.0078 



106.22 

5.26 

2 . 0552 

2 . 0557 

. 02 

200 . 01 

4.87 

2.0965 

2.0957 

- . 04 

297.05 

4.58 

2.1274 

2.1281 

. 03 

396.64 

4.394 

2.1491 

2.1502 

.05 

498.02 

4.267 

2.1649 

2.1661 

.05 


In summary, equation (18) agrees well with the observed changes in oscilla- 
tion amplitude for slow variations in tension and no damping. The energy losses 
seen in the retrieval runs are within that expected from the damping term in the 
retrieval tension control law. In the next section a more accurate expression 
is presented for calculating the energy losses from wire stretch damping. 




Figure 10 
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2.1.5 Rotational Energy Dissipated By Wire Damping Eorces - 

Rotation of the subsatellite will cause stretching of the wire except in 
the case of a perfect circular wobble. For an oscillation about an axis perpen- 
dicular to the wire, we can calculate the work done against wire stretch damping 
forces under the assumption that the center of mass remains fixed. The result 
of the derivation is 


W = b r 3 A 4 w 3 (23) 

L 8 32w J 

where W is the work done, b is the damping coefficient, A is the amplitude of 
the oscillation in radians, w is the frequency, and t is the time. The ampli- 
tude will, of course decrease with time so that the equation is valid only for 
small damping and short intervals of time. The average work W as a function of 
time is 


W = - b r 3 A 4 w 3 t (24) 

8 

Equation (24) would have to be used in a further integration to obtain the 
amplitude A as a function of time by setting the rate of loss of rotational 
energy equal to W/t. 

The details of the derivation of equation (23) are given in Section 2.0 of 
Technical Note TP86-002, January 1986 which is included as Appendix B of Quar- 
terly Report #6 for this contract. 

A simulation has been run using program ROTAT in order to verify equation 
(24) . In program ROTAT, the motion of the center of mass is not integrated so 
that the assumption used in the derivation is satisfied. The initial conditions 
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for the run are 6 = 2°, <p = ip = 0. The other parameters are moment of inertia I 
= 99 kg-m J , wire damping coefficient b = 326177 dynes/ (cm/sec) , wire tension F = 
3,990,530 dynes, and subsatellite radius r = 80 cm. The frequency w of the 
oscillation is 



.567862 rad/sec 


The amplitude A of the oscillation is .0349066 rad/sec. Putting all these pa- 
rameters into equation (24) gives 


W = 124.9 t 


Figure 11 shows a plot of the kinetic energy as a function of time during the 
12.5 second test run. Table 10 shows the kinetic energy at the two peaks. 


Table 10 


t 

(sec) 


Energy AE AT AE/AT 

(ergs) 


2.768548 194072.79 

8.301609 193386.87 -685.92 5.53306 -123.97 


This agrees with the calculated energy loss to an accuracy of .7 percent. 


Equation (24) can be used to make a more accurate estimate of the energy 
absorbed by the damping term in the retrieval control law. The damping term is 
given by equation (15) and the value of the coefficient is given in equation 
(16) . Using equation (17) with a 10 second period (w = .62832 rad/sec) , and 2° 
amplitude (A = .0349066) , the energy loss was estimated at .9 ergs/sec on the 
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basis of a crude analysis. Using equation (24) for the same parameters gives 
1.086 ergs/sec. In order to calculate the energy absorbed during the retrieval 
it would be necessary to integrate equation (24) over the whole retrieval. This 
has not been done for any of the cases. The assumption that the center of mass 
does not move would contribute some degree of error to the calculation. The 
major energy loss during retrieval is described by equation (18) . Equation (24) 
could be used to estimate the effectiveness of a reel control algorithm used to 
damp rotations of the subsatellite. 


2.1.6 Attitude Damping Using Active Tension Control - 

Extensive analysis has shown that the z-axis attitude control thruster on 
the tethered subsatellite cannot control small amplitude oscillations about the 
x and y axes. Since the coupling with these motions increases with amplitude, 
some control can be achieved at large amplitudes. Wire stretch damping is more 
directly coupled to the rotation but is also ineffective against small amplitude 
oscillations. It is possible in principle to use either tension control or the 
in-line thruster to control the torque on the subsatellite in a way which is 
effective against small amplitude oscillations. Both techniques require active 
control based on knowledge of the attitude of the subsatellite. Tension control 
using the reel motor on board the Shuttle would have to take into account the 
propagation delay of' tension signals along the wire in order to insure that the 
control is properly phased with the rotation. With the in-line thruster there 
is also a phasing problem since the subsatellite must move enough to stretch the 
wire and alter the torque applied by the wire. The thruster itself applies no 
torque since the line of action is through the center of mass of the subsatel- 
lite. For short tether lengths either method should work well. 


For thruster 
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control the delay time depends on the longitudinal oscillation period of the 
subsatellite at the end of the tether. There would be some immediate change in 
tension as a result of the inertial mass of the tether. This effect depends on 
the mass per unit length and propagation velocity for tension waves along the 
tether . 

A detailed analysis of the technique of attitude damping using tension 
control is presented in Section 3.0 of Technical Note TP96-002, which is in- 
cluded as Appendix B of Quarterly Report #6 for this contract. 


2.1.7 Attitude Damping Using A Lever Arm - 

Figure 12 shows a subsatellite with a lever arm between the wire and the 
subsatellite. In this analysis the lever arm is assumed to be long enough that 
it remains aligned with the wire (Figure 12a) . The lever arm could be attached 
to the subsatellite by a ball and socket joint or a universal joint having two 
perpendicular rotation axes. It is assumed that the joint is capable of apply- 
ing a damping tqrque that is proportional to the angular velocity of rotation of 
the joint. This arrangement is capable of damping any oscillations about the x 
and y axes of the subsatellite. It is ineffective against components of angular 
velocity that are parallel to the wire. This type of motion which is a spin of 
the subsatellite can be controlled by the z-axis attitude thruster. 
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This system is capable of damping a circular wobble which cannot be damped 
by any of the other damping systems described in this report. A circular wobble 
consists of oscillations about the x and y axes of the subsatellite that are 90° 
out of phase with each other. The component of angular velocity along the body 
z-axis is assumed to be zero since it can be eliminated by the z' -axis attitude 
control thrusters. Under this assumption, the damping joint must constantly 
rotate in order for the subsatellite to wobble regardless of the phase of the x 
and y components of the oscillation. 

For the sake of analysis, let us assume that the joint connecting the lever 
to the subsatellite consists of two axes, one of which is parallel to the body 
x-axis, and the other of which is parallel to the body y-axis. Let us assume 
that the wobble is described approximately by the equations 

0 K = A x sinw x t (25) 


and 


6 y = AyCOSWyt (26) 

where 0 X is the rotation angle about the x-axis and 0 y is the angle about the y- 
axis . For small amplitude oscillations equations (25) and (26) are a good de- 
scription of the motion. The frequencies of the oscillations are: 

Wx - vWlx ( 27 ) 


and 
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Wy = y/rE/Iy 


(28) 


where r is the distance from the center of mass to the attachment point, E is 
the wire tension, and I x and I y are the moments of inertia about the x and y 
axes respectively. The critical damping coefficients are: 


b x — 2 1 x oj x 


(29) 


and 


by = 2 IyOJy (30) 

I f I x = Iy and A x = Ay then the motion is a circular wobble. If the moments of 
inertia are unequal, the motion will alternate between circular or oval and 
linear wobble. The damping joints are effective against the x and y components 
individually and there is no requirement that the angle 6 between r and F be 
changing as long as the motion is not a pure spin. 

The lever arm may remain aligned with the wire either as a result of being 
long or as a result of the damping coefficient being low enough that the torque 
exerted on the lever by the subsatellite is small. The lever does not need to 
remain aligned with the wire for the system to provide damping. The configura- 
tion shown in Figure 12b will work also, but is somewhat more complicated to 
analyze. The only requirement for damping is that the joints rotate and thereby 
dissipate energy. 

The effect of a lever arm with damping joints is very similar to that of 
attitude damping thrusters about the body x and y axes. Some simulations of the 
effect of x and y attitude thrusters are given in Section 2.3.10 of Quarterly 
Report #5 , December , 1985 . 
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2.1.8 Comparison Of MMC And SAO Simulations Of Retrieval With Rotation Of The 
Subsatellite - 

As a further test of equation (18) comparison runs have been done using the 
simulation programs at MMC and SAO. The two programs were developed indepen- 
dently and each contains features not available in the other. Both model rota- 
tion of the subsatellite and retrieval with a tension control algorithm. The 
MMC program contains a model of the response of the reel control mechanism which 
is not available in the SAO version. In order to establish common ground for 
the comparison it was agreed to use as simple a case as possible that would show 
retrieval with rotation of the subsatellite. More complicated comparisons could 
be done later if needed once agreement was obtained on a basic reference case. 
Two or three runs were done with each program before establishing a common case 
for comparison. The final result was a retrieval run which treates the Shuttle 
as a point mass, and includes no model of the reel mechanism, a one dimensional 
rotation of the subsatellite, and a retrieval tension control law containing 
three terms involving damping, gravity gradient compensation, and a prescribed 
length vs. time profile. 

The runs done at SAO, in addition to being useful for comparison, are also 
useful for completing the analysis of an earlier case that needed to be rerun 
because of a program bug. Three simulations have been done at SAO for the 
comparison of the MMC and SAO analysis programs. The tension control law nor- 
mally used in the SAO simulation program is given by equation (2.1.16) on page 8 
of Quarterly Report No. 3, June 1985, for this contract. In the first simula- 
tion the program was modified to use the equation: 


f 


m*[3n J £ + 2fw c (£-£ c )] 


(31) 



Page 


65 


where : 

m* _ mim2 = reduced mass 

n = orbital angular velocity 

£ = distance between mass centers 

f = .9 = damping ratio 

w c = y / 3 fJ = control frequency 
£ c = - ■ ^ l 0 cos0 c sin0 c (0 C = 12°) 

The parameters for the run are : 


Tether length £ 0 
Shuttle mass mj. 
Subsatellite mass mj 
Shuttle altitude 
Subsatellite radius 
Moment of inertia 


20 km (deployed upward) 

101,614 kg 

550 kg 

296 km 

80 cm 

124 kg-m 2 


0 


The initial conditions for the run have the tether aligned with the local 
vertical and the subsatellite rotated 2° toward the forward direction (the di- 
rection of motion) . There is no initial angular velocity of the subsatellite 
(except for the orbital angular velocity) and no initial retrieval rate or 
tether libration velocity. 

The control law given by equation (31) differs from the control law cited 
in Quarterly Report No. 3 in that the factors accounting for deceleration of the 
subsatellite and the reduced gravity gradient for non-vertical tethers have been 
omitted. The effect of omitting these terms is to give a larger control ten- 
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sion , and more rapid retrieval at a larger retrieval angle. The control fre- 
quency w c for the damping term is taken as the in-plane frequency rather than 
the out-of-plane frequency as in Quarterly Report No. 3. 

A simulation has been done for 10,000 seconds of orbital time with output 
at 1 second intervals. The frequent sampling is necessary because the rota- 
tional period of the subsatellite is on the order of 10 seconds. The in-line 
thruster is never turned on in the SAO . Figure 13 shows the results of the 
first simulation. Part (a) is the tether length vs. time. The final tether 
length at 10,000 seconds is 117 meters. Part (b) shows the in-plane angle of 
the tether vs. time. The angle is larger than the commanded angle of 12° be- 
cause the tension is somewhat larger than required in this simplified control 
law. The effect of the subsatellite rotation can be seen as the tether length 
becomes shorter near the end of the run. Part (c) shows the tension vs. time 
and part (d) shows* the rotation angle of the subsatellite. The amplitude is 
almost 8° at the end because of the large tension ratio between the first and 
last points. In the current mission plan, the in-line thruster would be acti- 
vated to maintain adequate tension during the final stages of retrieval so that 
the growth of the oscillation amplitude would be halted before it reaches the 
levels shown in Figure 13d. This run provides an opportunity to test equation 
(18) for the case of retrieval from a length of 20 km with a single-axis rota- 
tion of the subsatellite (no coning) . Table 11 shows the results of the compar- 
ison. The amplitudes of the oscillation shown in the second column are obtained 
by quadratic interpolation between the output points at 1 second intervals, and 
the tension values are crudely interpolated between output points. The ampli- 
tudes in the fourth column calculated using equation (18) show good agreement at 
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the beginning of the run with the percentage error increasing with time. The 
energy losses have been shown to be within the limits expected from the damping 
in the control law. 


Table 11 


t 

e 

Tension 

d 0 (T„/T) v* 

%Error 

£ 

(sec) 

(deg) 

(dynes x 10 6 ) 

(deg) 


(km) 

0 

2.0000 

5.79590 



20 . 000 

5.1 

2 . 0016 

5.77063 

2.0022 

- . 03 

19.999 

105.0 

2 .0428 

5.3209 

2.0432 

- . 02 

19.874 

203.2 

2 . 0756 

4.9779 

2.0775 

- . 09 

19.587 

402.3 

2.1234 

4.5366 

2.1263 

- .14 

18.752 

1003.9 

2.1692 

4.1563 

2.1734 

- .19 

16.439 

2000.5 

2.4665 

2.4633 

2.4770 

- .43 

12 . 243 

3009.4 

2.8691 

1.3327 

2.8882 

- .67 

7.241 

4009.7 

3.2370 

.8111 

3.2700 

-1.02 

4.320 

5014.5 

3.7649 

.4346 

3.8220 

-1.52 

2.484 

6015.4 

4.3465 

.2385 

4.4406 

-2.16 

1 . 376 

7010.3 

5.0073 

.1309 

5.1591 

-3.03 

.763 

8006.9 

5.7997 

. 0695 

6.0439 

-4.21 

.415 

8990.0 

6.6805 

. 0371 

7.0708 

-5.84 

.224 

9961.0 

7.6694 

. 0197 

8.2831 

-8.00 

.120 


The first run, although useful 

as a 

test 

of 

equation (18) cannot 

be com- 

pared to any of 

the MMC runs since the 

MMC 

tension control law contains an 

additional term. 

In order to make 

sure 

that 

the 

retrieval profile is 

correct 

before doing any 

further rotational 

dynamics runs 

(which consume a lot 

of com- 


puter time) , a test has been done for retrieval without rotation modelled. The 


tension control law used for this run is: 


f c = m* [ 3f2 2 £ + 2fw c (£ — £ c ) + w c 3 (£-£ c )] 


(32) 


where : 



io 

lo 

a 

and the other parameters are the same as in equation (31). The orbital parame- 
ters and initial conditions are the same as before except that the subsatellite 
is treated as a point mass. Since there is no rotational dynamics, the integra- 
tion proceeds rapidly because there are no other short period effects to be 
followed. The integration has been run for 10,000 seconds of orbital time with 
output at 100 second intervals. Figure 14 shows the results of the simulation. 
Part (a) is the length vs. time. The final length is 446.6 m at 10,000 seconds. 
The tether stops retrieving at around 5600 seconds at a length of 2.227 km and 
redeploys slightly to a length of about 2.257 km at around 6000 seconds before 
resuming the retrieval . This is presumably because the retrieval was running 
ahead of the commanded profile specified by the variable £ c . Figure 14b shows 
the in-plane angle vs. time. The behavior here differs markedly from that seen 
in Figure 13b. The initial overshoot to about 30° is counteracted by the slow 
down in the retrieval rate and the tether swings to about 10° in the opposite 
direction. Figure 14c shows the behavior of the tension vs. time. The tension 
does not decrease monotonical ly with time. 

The results in Figure 14 look closer to the results shown in the first MMC 
simulation. However, differences remained. The MMC simulation contained a 
model of the response of the reel mechanism which is not included in the SAO 
model. When the reel effects were removed from the MMC simulation, agreement 
was achieved between the shapes of the curves from the SAO and MMC simulation 
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programs . 
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Having achieved satisfactory agreement on the control laws the retrieval 
run has been repeated with the rotation of the subsatellite included. Figure 15 
shows the results. Part (a) is the length vs. time. The final tether length is 
445.8 meters (446.6 m between mass centers with r = 80 cm) . Part (b) is the in- 
plane angle vs. time, and part (c) is the tension. As expected, Figures 15a, b, 
and c, are virtually identical to Figures 14a, b, and c. Figure 15d is the 
rotation angle vs. time. The shapes of the curve seems to agree with that of 
the MMC simulation at least up to the point where the in-line thruster comes on 
it the MMC simulation. Figure 15e shows the rotation angle during the first 500 
seconds, and Figure 15f shows the last 500 seconds. The difference in frequency 
results from the decrease in tension during the run. Figure 15g shows the 
rotational kinetic energy as a function of time. The difference in amplitude on 
alternate peaks results from the fact that the rotational kinetic energy is 
measured in inertial space. The angular velocity of the subsatellite adds to 
the orbital angular velocity on the forward swing and subtracts from it on the 
backward swing. 

The rotation amplitude has been analyzed using equation (18) the results 
are shown in Table 12. The agreement is good at the beginning, and the energy 
losses are within what can be expected from the damping in the tension control 


law . 
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Table 12 


t 

6 

Tension 

0 o (T o /T) 1/4 

%Error 

£ 

(sec) 

(deg) 

(dynes x 10 s ) 

(deg) 


(km) 

0 

2.0000 

5.79590 




5.14 

2.0011 

5.776 

2 . 0017 

- .03 

19.999 

104.62 

2.0319 

5.395 

2.0362 

-.21 

19.872 

201.97 

2.0669 

5.074 

2 . 0676 

- .04 

19.572 

404.96 

2.1175 

4.599 

2.1191 

- .07 

18.641 

1003.84 

2.1646 

4.191 

2.1689 

- . 20 

15.971 

2006.70 

2.5593 

2.126 

2.5699 

- .41 

. 10.950 

3004.40 

3.2368 

.822 

3.2591 

- .69 

5.880 

4010.40 

3.8100 

.4190 

3.8571 

-1.24 

3.319 

5017.58 

3.7785 

.4247 

3.8441 

-1.73 

2.357 

6011.75 

3.3976 

.6358 

3.4752 

-2.28 

2.257 

7001.66 

3.7532 

.4187 

3.8578 

-2.79 

1.845 

8003.43 

4.3739 

. 2213 

4.5244 

-3.44 

2.268 

9000.72 

5.1518 

.1111 

5.3750 

-4.33 

.697 

9995.02 

5.5797 

.0769 

5.8929 

-5.61 

.447 
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2.2 Slack Tether Studies 


, During the reporting period SAO has studied the situation in which the 
deployer reel jams during the retrieval maneuver, leading to loss of tension 
with the tether and satellite moving toward the Shuttle. Specifically, SAO has: 

Studied, in considerable detail, the dynamics of the tether following reel 
jam up to loss of tension; the tether velocity on loss of tension is 
derived. 

Implemented the results of this study as initial conditions for the fully 
slack tether simulator SLACK3 , and made sample runs. 

Calculated the closest approach of the satellite to the Shuttle on a free- 
orbit assumption, and hence the maximum amount of tether available to 
impact on the Shuttle. 1 


More detailed study of specific cases is planned for the next reporting 
period. 


2.2.1 Reel Jam: Dynamics Up To Loss Of Tension - 

The dynamics of the system in the reel jam case are trivially different 
from the model in the simulation program SLACK3 -- the only major difference is 
that the free end tether segment is replaced by a massive satellite. However, 
SLACK3 (as detailed in previous reports) of necessity deals only with fully 
slack tethers. It cannot follow the process from the initial impulse (tether 
break or reel jam) up to the time tension is lost, since elastic forces are 
dealt with only indirectly. Thus, we need to determine separately what happens 
during the initial loss of tension process, so that we may prepare "initial 
conditions" to start the simulation in SLACK3 . 

In this section we apply general wave propagation techniques to the problem 
of the initial dynamics following reel jam, and derive an explicit algorithm for 
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computing the tether velocity upon loss of tension. Generally (though see Sec- 
tion 2.2.3) , the reel jam will generate an impulsive wave (step function profile 
of strain) which propagates up and down the tether decreasing strain further at 
each pass and reflecting off the end masses until tension is lost. 

In Section 2.2.2 we apply separation of variables to assist in the above 
examination; unfortunately separation of variables breaks down at the moment of 
loss of tension, and we must combine the two approaches. In Section 2.2.3 we 
make an approximate calculation which shows that for the given reel control law 
and typical system parameters, we will be in the simplest regime, where the 
tether loses tension on the first pass of the impulse, and hence we need not be 
too concerned with the fact that we have greatly simplified the end reflections. 


2. 2. 1.1 Pre-Jam Assumptions - 

To facilitate obtaining a solution we make the following explicit assump- 
tions : 

- Uniform Tether Velocity V 

- Uniform Tether Strain e 

- At t = tj = 0, the reel/Shuttle turns into an infinite mass at rest, 
and the satellite becomes an infinite mass moving at velocity V. 

The tether velocity V, due to the reel motion under the retrieval control law, 

we also call the "reel velocity". 

These assumptions ignore several potentially important factors: 

- V, e variations due to gravity gradient forces on a finite mass tether 

- V, e variations due to reel acceleration 

- structure of attachment at Shuttle 

The V and e variations probably do not strongly effect the overall result we 
shall obtain. Eor short tethers, one or even twenty kilometers, the satellite 
mass dominates the tether mass, and the gravity gradient induced variations will 
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be small. The time scale of the retrieval law is long enough (hours) compared 
to the tether dynamical time scales that any variations in V or e due to reel 
motion will rapid equalize along the tether. The attachment at the Shuttle is 
more problematical; any reflections will clearly not be clean. Our simplifica- 
tion of the attachment is partly due to practicality: the structure is not 
adequately defined, and its inclusion would greatly complicate the analysis. 
Also, we show in Section 2.2.3 that the retrieval law is such that we will 
almost always be in a regime where no reflections of loading waves off the ends 
will be expected. 


2. 2. 1.2 Problem Parameters And Variables - 


Throughout 
V, e 
EA, H 
E 1 A (= Cv) 

L 

m 

s 

n 


the analysis we shall use the following notation: 

: initial (uniform) velocity, strain 

: tether properties, typically 10 s kg m/s 2 , 10" 2 kg/m 

: tether damping; use = 0, actually ~ 200 kg m/s 

: tether natural length deployed at tj = 0 

: subsatellite mass 

: tether natural coordinate 

: orbital angular velocity, typically 10~ 3 s' 1 
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2. 2. 1.3 Velocity Scales - 

We may define several velocities in the problem, all of which will be of 
importance : 

c : speed of sound in the tether; speed of wave propagation 

V : initial velocity (reel velocity) 

Vrec : Recoil velocity of cut tether; a function of e. 

We can estimate the range of these velocities, information which will be 
useful later. The simplest is the sound speed c, which depends only on the 

tether mass properties : 

c = \/EA/p 

Typically, c ~ 3 x 10 3 m/s. 

The recoil velocity is simply V REC = ec. To estimate this we need an esti- 
mate of e at the moment of reel jam. From the retrieval law in the next section 

with 9 c; 10°, T O' 2.7 Ci 3 mL . Take O ~ 10' 3 s' 1 , m ~ 50 0 kg, and L ~ 10 3 m. 

Then T ~ 1.35 kg m/s 3 (Newtons) and hence e = T/EA ~ 10' 5 . Then the typical c 

~ 3 x 10 3 m/sec implies V REC ~ 3 x 10' 3 m/sec. Similarly if we take L ~ 2 x 

10 4 m (20 km) we get T ~ 39 kg m/s 3 , e ~ 4 x 10' 4 , and V RBC ~ 1.2 m/sec. 

The control law gives a reel velocity depending on the retrieval angle and 
the current length, V = (3/4) Lflsin (29 ) . With 6 = 13°, V tv 0.3 L (km) . 

In summary: 


L = 1 km, 

Vr E c ~ 0.03 m/sec 

V ~ 

0 . 3 m/sec 

L — 2 0 km , 

Vr E c ~ 1 m/sec 

V ~ 

6 m/sec 


Note that the reel velocity V is greater than the recoil velocity V REC by about 
an order of magnitude, and the sound speed c is some three orders of magnitude 
higher still. 
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2. 2. 1.4 Retrieval Law - 

The control law assumed for retrieval is 

L = -a L 

where 

a = ^ Cl sin (2 6) 

The equilibrium behavior generated by this law is retrieval at a constant angle 
9 ahead of the vertical (for an upward deployed tether) . 


The tether tension is given (essentially) by two factors, (1) gravity gra- 
dient force on the satellite, and (2) reel acceleration/deceleration by the 
control law. We have 


T gg = (3 m tl J co s 0 ) L 

and 

d-d 

Tcl = -m L = -m — (L) = -m — (-qL) = a m L = -a 3 mL. 

Note that these have the same variation with L and Cl. The total tension is then 


T = [3O 2 cos0 - — n 3 sin 3 (2 0)] mL — 3 n 3 mL [cosO - 3sin 2 0cos 2 0] 
16 


For a typical value of 6, 13°, 

Tgg/TcL 


3cos0 

9/16 sin 3 (20) 


25 


In this case, the gravity gradient induced tension strongly dominates. However, 
if we take 9 = 45°, then T 99 /Tcl is only 3.77. 
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2 . 2 . 1 . 5 Gedanken Experiment V = V rec - 

Energy considerations can specify the magnitude of the tether velocity 
change (by conversion of elastic into kinetic energy) but do not tell us the 
direction of the increment. This simple energy balance also does not apply 
after the loading wave reflects from the ends since energy will then be trans- 
ferred to the Shuttle or satellite. To resolve these ambiguities and to pave 
the way for further analysis, we perform a simple thought experiment: 


First, transform the reference frame so that tether is at rest. 

s HUTTLC- 

V = o 
STKai*' £. 

t 

h»v 

We refer to this frame of reference as the "Tether Frame", and the frame in 
which the Shuttle is at rest as the "Shuttle Frame." Now, cut the tether at the 
moment the reel jams. 

CUT TIT HER 

h* C 

L 0 fc l rl <? vv /> '/£• 

f — > v 

The free tether is now the familiar cut-end case. A wave of loading will 
propagate down the tether at velocity A = c, and the tether in its wake will be 
moving to the right at velocity 

Vrec = e c 

If the Shuttle is moving at precisely this speed, V = V BEC , it makes no differ- 
ence if we have cut the tether or not; i.e., we have found the solution for 






V =o 


■Q 


s*TtL l ire 
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reel- jam in the case V = Vrec. 

In this case , V = Vrec, the tether will relax to e = 0 and be hanging dead 
in space — except that the satellite m is still coming in at V and will collect 
tether as it comes . 

This special case, V = Vre C , divides reel jam cases into two regimes, V > 
Vrec and V < Vre C . The two cases each have unique features and will be treated 
separately. 

2 . 2 . 1 . 6 Fast Retrieval , V > V rec - 

The only difference from the V = Vrec case is that the Shuttle is moving 
faster than the tether which has gone slack and catches up with a portion of it 
before the slackening wave reaches the satellite. We are left, at t = tj + L/c 
with: 



or, transferring to the original reference frame: 
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Note that there might be some change in the velocity of the mass m during the 
time L/c that the wave takes to release tension, but: 

- This time is short. 

- We are supposing that the tether tension basically balances the gravity 
gradient force. 

- If we worry about this then we should also take into account Coriolis and 
drag. 


2 . 2 . 1 . 7 Slow Retrieval , V < V rec - 


In this case, our cut tether analogy soon leads to: 



I.e., the tether separates from the Shuttle. Of course, in the actual reel- jam 
the tether remains fixed to the Shuttle which moves inexorably at V. We will 
have 

o Ak It' 6 w « V|f 
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2 . 2 . 1 . 7 . 1 Jump Conditions - 

In Gullahorn and Hohlfeld (1986) we derived jump conditions across the 
boundary between slack and taut regions. We now need the analogous conditions 
when the tension is positive in both regions. Following Appendix J of Gullahorn 
and Hohlfeld: 

Mass Conservation : No forces are involved in this derivation, so it is un- 

changed from slack case: 


A*i (A - v^ = p 2 (A - v 2 ) 

Momentum Conservation : The derivation in Gullahorn and Hohlfeld is unchanged up 

to : 

-Po a (Vj-Vi) = T 2 -T x 


Here the p's are linear densities, p 0 = Pi(l+€i) = p 3 (l+e 2 ) is the "natural" 
density. A is the physical velocity of the boundary; A = (A - v x )/(l+€i) = (A- 
V 2 )/(l + e 2 ), the velocity in "natural" coordinates. Using 
Vi = 0 v 2 = V 


ei - £ i — Cf 

we get, with Hooke's law (EAe = T) , 
pi A = p £ (A-V) 

-Pi A V = EA (£ f -ei) 
The first of these equations gives: 


leading to 


(Pi-p { )A = -p f V 


A 


Mt v _ A/ (1 + gf) 

Pt-Pi A/ (1 + e f) ~ 1/ (1 + £ i) 


= 1 + €l v 


Cl — £ f 
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Then the second equation becomes 



Noting that the tether undergoes compression, e f < ei, and that V > 0, 

€1 - € f = V/C 
or 



i V 1 


1 V ] 

ft = £i 

£ i c 

= « l 

V RE c_ 


Also note 


or 


A = 1 + e -~ • c(£i-e f ) = (l + fi)c ft* c 

£i - £f 


A = I" Vi ^ c 
1+ei 

I.e., the boundary propagates as expected, at the speed of sound. 


In summary, up to t + L/c we have 



As a check, note that we can obtain the same result for Cf from a kinematic 
argument. If Xi = €iL and X* = CfL are the physical lengths of the tether before 
and after the loading wave passes (L is the natural length) then at t = L/c, 
when the wave has just reached the satellite, we know the distance between the 


Shuttle and satellite has closed by Vt , giving: 
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Xf = Xi — V — . 

c 

=> £ t L = Cj. L - - V 

c 



as above . 

Note that 

- The tether is still taut. 

- Therefore, the stress wave can reflect from the subsatellite. 

- The Shuttle and satellite are still closing at velocity v = V. 

These imply: 

A. We need to look at reflection of stress waves. Note that (assuming 
clean reflections and no damping) the stress profile will always be a 
step function with one step . 

B. The Shuttle and satellite will eventually close to a distance allowing 
complete slackness at t a = (c L)/V. The time to propagate a wave along 
the tether, i.e., the time per bounce, is St = L/c . Thence, the number 
of bounces before the tether goes slack is about 

N = — = <L / V = ££ _ v rec 
8 St L/c V V 

Of course, pending analysis: 

- there may be slack regions before t„ 

- even after t a there is no guarantee that one will achieve uniform 
slackness. 
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2. 2. 1.7. 2 Impulse Reflection And Evolution Past T = L/C - 

Attempts to apply general reflection considerations lead to some confusing 
subtleties . 

Instead, we have solved explicitly a sample problem using separation of 
variables, and evaluated the (series) solution with a simple computer program. 
The derivation is below in Section 2.2.2. One can scale the problem so that the 
only free parameter is V/V BEC (or simply V, since V BEC is scaled to 1) . 

Sample results for V/V REC = -0.3 are on the next page. In general, we may 
induce 

1. The wave front propagates at c, reflecting instantaneously and continuing 

at c . 

2. The wave is always compressive by the amount 6e = V/C (= £i) . I.e., 

nothing unexpected happens at the reflections. 

3. The velocity is 

v = 0 to the left of wavefront 
V to the right of wavefront 

Alternately, when the wavefront passes a mass element there is an increment 

in its velocity of magnitude 

6v = | V [ 

The increment is in the direction the wavefront travels. 

Everything is clear-cut up to loss-o f- tension . Now we must consider the 
final pass of the wave through the tether, during which tension is lost. 

F irst t consider a case where the wave is about to reflect from the left 
(fixed) end, and would thereafter have negative tension: 












After reflection we will have a situation similar to the V < Vd^ treated 
above : 




The velocity jump, as before, will be "Vrec," but in this case the V REC appro- 
priate to the remnant strain , 

Vf — -V + "Vrec" = -V + 6 fC 

y 

Noting that e t < 6c by definition, and that 6c ** — (= initial V„ EC ) , we confirm 

c 

that 

v t < — V + — c = 0. 
c 

All we can (generally) say is 

- V < v t < 0 

The above conclusions are for the situation with the moving boundary at the 
right: 
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and we get 



or, transferring to the desired Shuttle coordinates: 



v~o 


2. 2. 1.7. 3 Enumeration Of Final States For Decreasing Reel Velocities - 

As we decrease the initial (reel) velocity from some value larger than 
V RECl the loss-of-tension problem passes through distinct regimes as the number 
of reflections of the compressive wave before complete loss of tension in- 
creases. We shall examine explicitly the first few of these regimes to seek 
easily encapsulated regularity. Define a parameter 
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The initial configuration is: 



V* 


e ( 


F/A ££> 


v > 1 i.e. V > V REC . In this case the tether loses tension on the initial 
pass of the compressive wave. No reflections need be considered. . 


V f = V - ec = V - V REC 

= [l/ - 1] £ C 


= [‘ - ft 


An amount V»— of tether will be left crumpled on the Shuttle. 


— <«/<! The loading wave will reflect once, off the satellite. 


Here, as in all cases below, 


S € = — = !/£ 

C 


For the single reflection case, 


= € " 6e = [ x ~ 7^] e = cl_,/]€ 


Note that 0 < e t < . In particular, e t (i/=l)=0, and £ f (l/2) = e/2. The final 

velocity is: 


V f = V — v f = £ f c = [l—i/] ec 
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In this case, a portion of tether will be crumpled onto the satellite. 


i < * < i 

3 2 


There are two reflections, the final one off the Shuttle. 


V, = V — 


e { = e — 26 e = [1 — 2 v] e 

“ (t) - " (I) = I ' 

ffC = vec — [1 — 2i /] ec = [3v — 1] ec 


Note that Vf(i/ = 1/2) = ec/2, and V*(l/3) = 0. 


1 . . 1 

— < v < — 

4 3 


There are three reflections, the final one off the satellite. 
e t = [1~3i/]€ 


" (I) - • •«(;)■ 3 


V c = € t c = [1 — 3i/] ec 


Vf (i) - i ec V t 


©■ 


The trend is clear, and is plotted on the next page. 


Summary : The result is cleaner if we use — instead of v. 


l/i/ < 1 

1 < l/i/ < 2 

2 < 1/v < 3 

3 < 1/v < 4 


V f = (i/-l) ec 
V t = (l-i/) ec 
V t = (3i/-l) ec 
Vf = (l-3i/) ec 


C - 
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Figure 17. The final tether velocity V f as a function of the reel velocity at 
the moment of jam, V. When scaled by the elastic recoil velocity Vrec = £c (the 
recoil velocity of a cut tether), this simple functional relation results. This 
result ignores gravity gradient effects (small during the few seconds of re- 
coil) ; any initial tension variations due, for instance, to variations in reel 
velocity prior to jam; and assumes an undamped Hooke's law tether. 
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This generalizes to: 


2n < - < 2n+l V ( = 

V 


2n+l < - < 2n+2 V f = 

V 



| (2n+l) i/-l| 
|l- (2n+l) 



£C 

£C 

ec 

2n+2 


2. 2. 1.8 An Algorithm For The Initial Velocity - 

From the above discussion we can now write down an explicit algorithm for 
computing the final tether velocity upon loss of tension after the loading wave 
has finished propagating up and down the tether reflecting at either end. As 
input data we need only the initial reel velocity, the initial strain, and the 
tether properties . 

1. ' Compute v V/ec, where: 

V = initial velocity, at reel-jam 
e = initial strain , e = T/EA 
c = speed of sound, c 2 = EA//i 

2. Compute n = w hei"® [•] is the greatest integer function. 

3. IF 2n < — < 2n+l THEN 

V 

Vj = j(2n + l)i'-lj ec 


ELSE 


V f = |l- (2n+l) i/J ec 
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2.2.2 Separation Of Variables Applied To Reel Jam - 

We apply standard separation of variable techniques to provide an explicit 
series solution of the reel jam problem. This solution applies only up to the 
first loss of tension since it relies on the wave equation being satisfied in 
the body of the tether; as soon as tension is lost, the wave equation does not 
apply uniformly. It appears that the series could be evaluated in closed form, 
but to reduce both effort and possibility of error we have chosen simply to 
evaluate it directly. 

We consider the problem in the coordinate system comoving with the tether 
and satellite at the moment of the jam (t=0) . The physical coordinate is x, and 
the natural coordinate along the tether is s. The satellite is fixed at x = s = 

0 . (Note that we are reversing the ends from the convention used in Section 
2.2.1.) The Shuttle, with s = L, is initially at x = (1 +£i)L, where L is the 
tether deployed natural length at the moment of reel jam and £* is the initial 
strain. We ultimately seek the strain and velocity profiles £(s,t) and v(s,t), 
but we shall work with the displacement u(s,t) = x(s,t) - s. Then the tether 
wave equation (EQ) , initial conditions (IC) and boundary conditions (BC) are: 

EQ u tt — c 2 u„ = 0 0<s<L,t>0 

IC u(s,0) = £i s 0<s<L 

u t (s , 0) = 0 

BC u (0 , t) = 0 t > 0 

u(L,t) = € t L - Vt = q (t) , say 

The first step in solving these equations is to make the BC homogeneous, 

1. e. with zero right hand side. Introduce 

4 >{ s,t) = (s/L) q (t) 
and then transform to a new dependent variable u, 
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u(s, t) = u(s, t) 

- 0(s,t) 

The 

system then becomes : 


EQ 

o 

II 

m 

a 

r« 

0 

1 

V 

0 < s < L 

IC 

u (s , 0) = 0 

0 < s < L 


u t (s,0) = (V/L) s 


BC 

u (0 , t) = 0 

t > 0 


u (L , t) = 0 



Note that in general EQ would now be inhomogeneous, with right hand side 
-(s/L)q"(t) , which would lead to further complications. We have this simple 
system because the 'forcing' motion of the end point (Shuttle) is at constant 
velocity . 

The above is now a simple textbook exercise, and separation of variables 
gives a solution of the form: 




n=l 


An sin^ ^ + B„ cos( 


njrct 



This formula incorporates the BC. We must use the IC to determine the 
coefficients A„ and B„, completing the solution. 


oo 

The first IC gives 0 = u(s,0) = B„ sin (n7rs/L) . This can only be true 

n=l 

for all s € [0,L] if B„ = 0 for each n. 

OO 

The second IC gives (V/L) s = u t (s,0) = A„ (n7rc/L) sin (n^rs/L) , 0 < s < 

n=l 

L. The standard trick of multiplying by sin (m7rs/L) and integrating from 0 to L 
gives, after some algebra, A,, = 2 (VL/c) (-l) m+1 / (m?r) J . 

The solution is then: 


2VL 1 (~l) n41 



nffct 

L 



u(s, t) 
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and the physical displacement 

u (s | t) = u(s,t) + V>(s,t) 

, = u(s,t) + { £iS - (V/L) St >. 

We ultimately want the strain, £ = 3u/3s , and the physical velocity v(s,t) = 
3x/9t = 3(s+u)/9t = 3u/3t . These are 

/ 4 -\ V [ 2V 1 r-' (— l) n+1 . ( n?rct \ ( nffs \ 

"•■t> = 1 ' l 1 + [ 7 j ) co i — ) 

and 

, . . v 2 ' (-l)"* 1 ( nirct \ . ( nws \ ,,s 

v(s,t) = 2VL E- L -^r~ cos ( — J Sln ( — ) - 

We have written a simple program to evaluate the above expressions for 
e(s,t) and v(s,t). One inputs the problem parameters V, L, c, and some run 
parameters such as the number of terms to sum in the series. The output is 
examined (although it could be plotted) and with 100 terms the result is quite 
clear: a step function in € and v propagating at c, reflecting at the bound- 
aries. Sample results are shown above in Figure 16 and discussed in Section 

2. 2. 1.7. 

The problem can always be scaled to leave just one parameter. 

First scale s so that L=1 . 

Then scale t so that c=l . 

Then scale u so that £i=l. 

Noting that the second and third scalings imply V REC = 1, we are left with only 
one essential parameter that we can choose, the boundary (reel) velocity V/V R ec 
(= V under the above scaling) . 
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2.2.3 Implications Of The Control Law For Reel Jam Initial Conditions - 

It was seen in the above that, at least under reasonable simplifying condi- 
tions, the motion of the tether upon final loss of tension is dependent upon 
only one parameter, v = V/ec, where V is the tether's velocity at the moment of 
the jam (reel velocity) , € is the tether's strain at that moment, and c is the 
axial speed of sound in the tether, c 2 = EA//i. 

Given the particular form of the retrieval control law in Section 2. 2. 1.4, 
we can compute this factor i/ as a function of the system parameters (tether 
properties and satellite mass) and of one variable, the retrieval angle 6. In 
particular, v does not depend on the length at reel jam. For typical system 
parameters , v is comfortably greater than 1 for retrieval angles greater than 
about 2°; thus, in practical cases, we avoid the complexities of the poorly 
known reflection properties from the satellite and deployer. 


The reel velocity will be given directly by the control law. We also need 
the tether strain e . In computing this we ignore the reel acceleration and the 
tether mass. We may do the first because we have seen above that the gravity 
gradient tension dominates, strongly at moderate and small angles. It will turn 
out that it is at a quite small angle that the transition to the high velocity, 
v > 1 , regime takes place. The tether mass will be dominated by the satellite 
mass for the short tethers (say 1 km) of interest. 


The gravity gradient force is then simply 

T, g = (3 m n 2 cos ^ ) L 

leading to initial strain 

£ = Tgg/EA 


_ / 3 m fl 2 cos 6 N L 

_ V EA ) 


EA 
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The velocity is given directly by the control law: 

V = |L| = aL 

where 

a = (3/4) ft sin (20) 

Here, L is the length at reel jam, ft is the orbital angular velocity, and 8 is 
the retrieval angle. 

Combining the above expressions: 


^ ft sin ( 2 6 ) 



ec 

/ EA 

3 m ft 3 cos 8 

or 


V Z 3 

EA 


v = { } sinW 


Note that this does not depend on L. That is, as long as we are in steady state 
retrieval (that is the actual tether angle is equal to the control law angle) , 
the basic form of the loading after the reel jam does not depend on how close 
the satellite is at Jam. The actual velocity at loss of tension will vary with 
L, but this will be a simple scaling; since the distance to be covered by the 
recoiling tether is similarly scaled, the time scale for impact on the Shuttle 
is also independent of the length at jam. 

We can evaluate the above expression for v with typical parameter values ft « 
10' 3 s' 1 , m as 500 kg, p as 10" 3 kg/m, AE « 10 5 kg m s" 3 . This gives 

1 / ft) 30 sin (S) . 

Thus, for sin 9 > 0.03, or 0 > about 2°, v > 1 and we are in the "slack on 
first pass" region. That is, the tether will go slack during the initial propa- 
gation of the wave of compression without waiting for reflections from the sat- 
ellite and/or Shuttle and the complications they introduce. Thus, in practical 
cases, the end result does not depend on the simplifying reflection assumptions 
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made in Section 2. 2. 1.7. 2, and the analysis is more secure than if it did. 

2.2.4 Satellite Orbit After Reel Jam - 

At the moment of reel-jam, or more precisely at the moment the tether loses 
tension, the satellite becomes a freely orbiting object. In examining the early 
behavior of the system, that is up until the tether is brought into full tension 
again, this free orbit will provide a good approximation to the actual satellite 
motion. While the tether is slack at the satellite end, it will obviously not 
effect the satellite motion, and even if a portion of the tether is affecting 
the satellite the influence will not be great since the tether mass is so small; 
it will only be when the tether acts simply as a taut link with the relatively 
immovable Shuttle that the approximation breaks down. 

In particular, we may compute the closest approach of the satellite to the 
Shuttle in its free orbit. This will give us an estimate of the amount of 
tether (tether length - closest distance) available to entangle the Shuttle and 
deployment mechanisms. 

We define our standard corotating coordinate system centered on the Shuttle 
with the x axis radially outward and the y axis along-orbit. The retrieval law 
from Section 2. 2. 1.3 will then specify the initial values of the satellite posi- 
tion and velocity in terms of the length at reel jam and the retrieval angle: 
x = L cos0 

y = L sin 9 

x 


Y 


— (3/4)nsin (20)L cos# 

— (3/4) flsin (20) L sin0 
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The equations of free body motion are 
x = 3fl J x + 2 n y 

y = - 2 n x - D 

where the drag is represented by 

D = C D — p. v a 2 . 
m 

Here, C D is the coefficient of drag, typically taken as 2; A and m are the 
satellite area and mass; p» is the ambient atmospheric density.; and v a is the 
velocity through the atmosphere. 

Scale distance by the tether length L, £ = x/L and r) = y/L , and time by 
the orbital angular velocity, r = O t, and denote the r-derivative by () 
Then we obtain the following scaled equations: 

Initial conditions: 

£ = cos 8 

r) = sin 8 

= — (3/4) sin (20) cos 6 

rf' = — (3/4) sin {20) sin 0 

Equations of motion: 

V' = 3f + 2,' 
r>" = - 2 | ' - (D/L) 

where D = D / (l 2 . 

We see that there are only two fundamental parameters left: 

8 -- The retrieval angle. 

D/L -- Drag / Length at reel jam. 

Note that 6 only appears in the initial conditions , and D/L only appears in the 
equations of motion. Also, if drag is negligible, the length does not enter 
into the solution and the only available parameter is 8. 
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We can evaluate a typical drag for 295 km altitude. Using Co = 2 , p, = 3 x 
10' 14 g/cm 3 , O = 1.2 x 10' 3 sec" 1 * v, = RQ = 8 x 10 s cm/sec, m = 500 kg = 5 x 10® 
g, ‘A = n (80 cm) 3 = 2 x 10 4 cm 3 , we get D fti 0.01 km. Thus, even for L quite 
small, say 100 meters, D/L is small. 

We have written a brief program to solve the equations of motion using a 
simple Runge-Kutta integrator. One inputs the retrieval angle 9 and the 
drag/length parameter D/L. The program then tabulates (on the terminal) the 
position and velocity, and also the distance to the Shuttle, as a function of 
time. We have exercised the program and noted the closest approach (in dimen- 
sionless radius, p 3 = £ 3 + t} 3 ) for various 6 and D/L. 


9 

5° 

D/L: 0.01 

0.1 

0.5 

0.9972 

0.9972 

0.9971 

15° 

0.9744 

0.9741 

0.9726 

30° 

0.8983 

0.8950 

0.8766 

45° 

0.7699 

0.7493 

0.6089 

60° 

0.5691 

0.4675 

0.0580 


If we assume D = 0.01, then the D/L value tells us L. We can then scale the 
above dimensionless minimum approach distances to get actual closest approach 
and subtract from L to get the amount of tether available to entangle the de- 
ployer (given in meters) : 


6 

L: 

1 km 

100 m 

2 0 m 

5° 


3 

0.3 

0.1 

15° 


26 

3 

0.6 

30° 


102 . 

11 

3 

45° 


230 

25 

8 

60° 


421 

53 

19 


We should point out that not all of the above tabulated tether will actually 
impact on the Shuttle. These values simply represent the maximum amount that 
the satellite's orbit will allow in the region of the Shuttle under "optimal" 
circumstances. The results tabulated above encourage the use of smaller re- 
trieval angles. The time scales to closest approach are also of interest. With 
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D/L = 0.01 (i.e. for a 1 km. tether) these are, at 5°, 33 seconds; at 15°, 120 
s; at 30° , 270 s; at 45°, 575 s; at 60° , 1200 s. For larger D/L these 
numbers increase, but not very substantially. 

These results have also helped us in interpreting the results of simula- 
tions reported below. 


2.2.5 Implementation Of Reel Jam In SLACK3 - 

We have added an option to the SLACK3 program for running simulations of 
reel jam cases (as well as the previous cut tether cases) . The only modifica- 
tions required are in the routine SPECIFY which sets up the initial state, that 
is the state at the moment the tether completely loses tension (SLACK3, to allow 
high resolution models, assumes all tether segments are slack, with instantane- 
ous "hard" bounces when a segment goes taut) . Recalling that SLACK3 uses a menu 
system (with defaults) for selecting many parameters, the apparent differences 
to the user are : 

- The reel jam option must be chosen (the default run type is still a 
tether break case, but this could easily be changed) . 

- If the reel jam option is chosen, the program attempts to force the user 
to choose reasonable deployment angles, e.g. in the orbital plane and 
with the proper sign for retrieval (forward if deployed up, backward if 
deployed down) . Out-of-plane or wrong-sign deployment could be chosen 
by specifying the angles before specifying the reel jam option. 

- If the reel jam option is chosen, the program asks for the length at 
jam, rather than the cut length. 

The choice of other parameters and options remains basically unchanged; thus, 
one could, for instance, have a reel jam case with Shuttle rotation and acceler- 


ation . 
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The assumptions made and internal mechanics of the modification are: 

- We assume that the actual deployment angle is the angle in the retrieval 
control law of Section 2. 2. 1.4. That is, the retrieval maneuver has 
reached the steady state. Thus, reel velocity does not need to be input 
separately. (This could be changed if desired.) 

- The internal tether masses are each given a velocity computed with the 
algorithm of Section 2. 2. 1.8. 

•First, the retrieval angle and the tether length give the reel ve- . 
locity. See Section 2. 2.1.4. 

•The angle and length also give the tether tension (which is assumed 
to be the gravity gradient tension only, ignoring the reel accelera- 
tion) using the tension calculation already in the program. 

• Combined with tether properties, we now get the strain € and sound 
speed c, and hence the dimensionless parameter v . 

•Vf/ec, and hence v £ , are now calculated as in Section 2. 2. 1.8. 

- The final mass point, formerly the free end, 

•has the satellite mass added to it, 

•is given the original reel velocity. 

- As previously, we make the tether slightly slack to begin with, with 
some randomization of segment lengths and directions. This was origi- 
nally done to avoid sudden uniform tensionings and simultaneous bounce 
problems, and in the reel jam case obviates having to deal with the 
(very small) amount of tether already swept up by the Shuttle or satel- 
lite while the final slackening wave propagates along the tether. 


In addition to these changes specific to the reel jam option we have taken 
the opportunity to fix a bug in the boom motion routine (it could no longer deal 
with the case of na Shuttle thruster firings) , and to change some of the parame- 
ter defaults. The default altitude and ambient atmospheric density are now 295 
km and 2.7 x 10~ 14 g/cm 3 (appropriate to 295 km and 10 0 0°) . The default tether 
properties are now p = 0.1 g/cm, AE = 10 10 dynes, Cv (= AE ' ) = 2 x 10 7 dyne-sec 
(see Lorenzini, et al . , 1985, Section 2.3.2 for these values). 

Two sample runs were made, both with reel jam at one kilometer. They 
confirm, in general, the satellite orbit considerations of Section 2.2.4. 
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(Those calculations were actually made in response to the SLACK3 runs.) The 
retrieval angles were 13° and 45°. 

In the 13° case, the progress of the tether toward the Shuttle is only 
barely visible on the plots, and the satellite turns around and draws the tether 
taut at about 300 seconds. (As noted in previous reports this generally causes 
the program to fail.) Note that the time to closest approach at 15° given in 
Section 2.2.4 is 120 seconds; since the SLACK3 simulation, starts out with the 
tether slack by a few percent, it. is not surprising that the time to full 
tension is somewhat more than the expected value of twice the time to closest 
approach . 

The 45° case is more interesting. In Figure 18 we see a substantial amount 
of tether impacting on the Shuttle in the first 200 seconds or so, followed by 
avoidance as the infalling tether is carried along orbit by Coriolis force. The 
tether seems to avoid wrap-around of the Shuttle. Although at the initial 
tether speed of 0.8 m/sec some 160 meters could have impacted on the Shuttle in 
that time span, visual inspection indicates that a smaller portion of the avail- 
able 230 meters did so. Note that the total time up to re-tension is 1225 
seconds, about twice the 575 seconds to closest approach from Section 2.2.4. 

More information could be obtained from the 45° case by rerunning it with 
more tether segments, and by modifying the plot program to show only the near- 


Shuttle segments. 
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3.0 PROBLEMS ENCOUNTERED DURING REPORTING PERIOD 

The analysis of the string-like dynamics of the tether following a thruster 
activation was one of the items planned to be started during the current report- 
ing period. Because of a rearrangement of the priorities, the above mentioned 
item has been replaced by the comparison simulations between SAO and MMA com- 
puter codes during retrieval and by extra analysis of the damping of the satel- 
lite rotational oscillations by means of reel control. 


4.0 ACTIVITY PLANNED FOR THE NEXT REPORTING PERIOD 

In the next reporting period we will investigate the effect of the wire 
intrinsical damping on the satellite rotational dynamics. A reasonable model of 
the tether damping as a function of tether length will be implemented into our 
comptuer code and retrieval simulations will be performed. 

The simulation activity of the tether dynamics following a reel jam will 
also be continued. Simulations of some of the most critical situations, during 
the retrieval phase, will be run by making use of the SLACK2/3 computer code. 



